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We consider a one-dimensional lattice of Ising-type variables where the ferromagnetic exchange 
interaction J between neighboring sites is frustrated by a long-ranged anti- ferromagnetic interaction 
of strength g between the sites i and j, decaying as | i — j |~", with a > 1. For a smaller than a 
certain threshold ao, which is larger than 2 and depends on the ratio J/g, the ground state consists 
of an ordered sequence of segments with equal length and alternating magnetization. The width 
of the segments depends on both a and the ratio J/g. Our Monte Carlo study shows that the on- 
site magnetization vanishes at finite temperatures and finds no indication of any phase transition. 
Yet, the modulation present in the ground state is recovered at finite temperatures in the two-point 
correlation function, which oscillates in space with a characteristic spatial period: The latter depends 
on a and J/g and decreases smoothly from the ground-state value as the temperature is increased. 
Such an oscillation of the correlation function is exponentially damped over a characteristic spatial 
scale, the correlation length, which asymptotically diverges roughly as the inverse of the temperature 
as T = is approached. This suggests that the long-range interaction causes the Ising chain to fall 
into a universality class consistent with an underlying continuous symmetry. The e^^'^ temperature 
dependence of the correlation length and the uniform ferromagnetic ground state, characteristic of 
the g = Q discrete Ising symmetry, are recovered for a > ao. 
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I. INTRODUCTION 



The competition between a short-ranged interaction 
favoring local order and a long-range interaction frus- 
trating it on larger spatial scales is often used to explain 
pattern formation in chemistry, biology and physicsii^. 
The role of the long-range interaction is to avoid the 
global phase separation favored by the short-ranged in- 
teraction and promote a state of phase separation at 
mesoscopic or nano-scales. Thus, the long-range inter- 
action is not, in general, a small perturbation^i^i^i^i^ii^, 
but must be considered as precisely as possible. From a 
computational point of view, this means that the frus- 
trating interaction has to be accounted for by involving 
all the lattice sites in the computation, which in turn 
limits the actual system size that can be handled in 
e.g. Monte Carlo (MC) simulationsSiia^iiii^iHiii. Few 
exact results on multi-scale, multi-interaction^i^ systems 
are present - to our knowledge - in literature. For 
one-dimensional systems rigorous proof of absence of a 
phase transition in the pure long-range antiferromagnetic 
model has been obtainedi^. Besides, rigorous results con- 
cerning the ground-state phase diagram can be found in 
Ref. i3j. Regarding two-dimensional lattice models with 
restricted spin orientation and dipole-dipole interaction 
competing with ferromagnetic nearest-neighbor exchange 
interaction, Giuliani et al}^ showed that the ground state 
is periodic striped, while a zero-temperature reorienta- 
tion transition (from in-plane to out-of-plane magneti- 



zation) occurs at a given relative strength of the short- 
and long-range interaction when are both antiferromag- 
netic. Finally, a generalization of this periodic ground 
state in some continuum versions has been rigorously 
proved^^. 

In this paper, we perform MC simulations on a one- 
dimensional (Id) lattice with sites occupied by Ising-type 
classical variables assuming values aj = ±1. The nearest- 
neighbor sites interact by a short-ranged ferromagnetic 
interaction of strength J which favors the same sign for 
two adjacent variables (in the language of magnetism the 
exchange interaction favors parallel alignment of neigh- 
boring spins). In addition, any two variables located at 
sites i and j interact by means of a long-range inter- 
action of strength g decaying according to a power law 
I i — j |~" and favoring, instead, antiparallel alignment. 
In the present study, selected values of a > 1 and in the 
vicinity of 2 are investigated. This range turns out to 
be representative of the different physical regimes. We 
are aware of the apparently academic nature of j.) a one- 
dimensional model and of ii.) this choice of values for 
a. In fact, point charges interact via the Coulomb in- 
teraction, which has a = 1, while the dipolar interaction 
between two localized magnetic moments has a = 3. On 
the other side, imposing a mono-dimensional modulation 
to two- or three-dimensional arrangements of charges and 
spins (a symmetry often realized in experimentsii^i^SiiSi) 
produces an effective one-dimensional long-ranged inter- 
action potential with an effective value of a which can 
differ from 1 and 3 respectively. As an example, ele- 
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mentary magnetic moments arranged into stripes and lo- 
cated on a two-dimensional array of sites interact with an 
effective, one-dimensional dipolar long-range interaction 
which, asymptotically, is proportional to \ |~^~. Ac- 
cordingly, a systematic study for values of a in this range 
might reveal properties that can be used to explain physi- 
cally relevant situations, such as those represented by the 
two-dimensional system of stripes quoted above or sim- 
ilar models of frustration discussed in connection with 
electronic phase separation^. A Id model has great com- 
putational advantages compared to its 2d and 3d coun- 
terparts, such as the possibility of simulating lattices of 
larger linear dimensions, which in turn allows larger mod- 
ulation lengths than already reported^ '^^i^^'^^i"'^^'^"^ , which 
are, indeed, closer to experimental situations. Later in 
the paper, we will single out the relevance of our results 
for understanding realistic spin and charged systems. Be- 
sides, variations of the Id-Ising model including long- 
ranged potentials have been widely applied to biological 
problems22i, such as protein folding^"^ and helix-coil tran- 
sitions^. 

This paper is organized as follows: In Section we 
introduce the model and its known ground-state phase 
diagramSi and then present our main results on the oscil- 
latory character of the two-point correlation function, on 
the temperature dependence of the corresponding mod- 
ulation period and on the correlation length. These 
facts point to the persistence of the modulated struc- 
ture emerging in the ground state even if, strictly speak- 
ing, the on-site order is completely lost in the thermo- 
dynamic limit^^. In Section IIIIl we provide some argu- 
ments aiming at explaining, within an analytic approach, 
the cross-over from the g = Ising universality class to 
a continuous-symmetry behavior for g > and a < ao 
as well as the temperature dependence of some physi- 
cal observables in comparison with MC simulations. In 
Section llVi we provide a summary of the most relevant 
results and indicate possible directions for further work. 
Technical aspects of the MC simulations and of the ana- 
lytical computations are presented in Appendices. 

II. MONTE CARLO RESULTS 

A. The model and the ground state 

The Hamiltonian with Ising variables ai — ±1 on a Id 
lattice reads 

N 

H=-j5:.,.,,,+f 5:-^, (1) 

where N is the number of spins in the chain and {i ^ j} 
indicates a sum over all the couples in the chain; peri- 
odic boundary conditions (Ti+Ar = ai are assumed. The 
ground stated of this Hamiltonian is uniform for a > uq, 
where ao > 2 depends on the ratio J/g. For a < ao, 
the ground state consists of a regular sequence of groups 
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FIG. 1: (Color online) Ground-state phase diagram in the 
{oi,g/J) plane. For 1 < a < 2 the ground state always con- 
sists of domains (grey region). For q > 2, the crossover to a 
uniform ground state (white region) occurs when a exceeds 
the threshold value qq indicated by the solid line. MC cal- 
culations have been performed for values of a between 1.6 
and 3.2 and g/ J = 0.4 (horizontal red line). Inset: zoom of 
the region where MC simulations have been performed in the 
{a , J/g) plane, J/g — 2.5 (horizontal red line). 

of h adjacent spins with positive {(Jj = +1) and nega- 
tive {aj = —1) orientation. The zero-temperature phase 
diagram in the {a,g/J) plane is schematically reported 
in Fig. [1] The inset zooms into the region of the pa- 
rameter space {{a, J/g) plane in this case) in which 
MC simulations have been performed: J/g = 2.5 and 
q; = 1.6...3.2. The main thermodynamic observable we 
address is the two-point correlation function at temper- 
ature T and fixed a: 

Ca(r) = {{<7j+ra,)j)T (2) 

and its Fourier transform iS^ (q) (commonly named struc- 
ture factor): 

-t-oo 

E ((^.+.^.).)Te-^^^ (3) 

r— — oo 

As the system cannot be assumed to have translational 
invariance, an average over the lattice sites j is needed 
in ([21) and henceforth); (. . . )t denotes the 
thermal average. The physical quantities computed with 
the MC approach actually correspond to the double av- 
erage (...) = ((...)i)T. 

The lowest-energy spin profiles are known to be square 
waves Sq(fcoj), with a modulation period 2h — 27r//c(j^. 
The total energy can be parameterized with h by in- 
serting the square profile into the Hamiltonian ([1]). The 
ground state equilibrium value of /i ~ let us call it hgs, 
corresponding to kgs - is then determined by minimizing 
the resulting energy (jB3p with respect to h. hgs depends 
on a and J/g: some values are reported in Fig. O The 
two-point correlation function for a generic square-wave 
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FIG. 2: (a) Correlation function for a = 1.8 (domain ground 
state) at different temperatures for L = 1000 and J = 2.5. 
(b) Correlation function for a — 3.0 (uniform ground state) 
at different temperatures for L = 1000 and J — 2.5. Inset: 
Reminiscence of the competing dipolar interaction (see text). 

profile reads (see Appendix |B] for details): 
1 ^ 

{crj+rO-j)j = — ^Sq(fco(j + r))Sq(A:oj) 
1 °° 

= 2 X! "™ ^ Tr(/sor) , 

m— 

(4) 

where ko — ir/h, km = (2m+l)-fco, a,„ = 4 [7r(2m + 1)] ^ 
and Tr(fcor) is a symmetric triangular wave of period 
2h. According to dU evaluated in h = hgs, the 
ground-state structure factor only takes non-zero values 
in the points located at q = ±(2m -I- 1) ■ kgs for which 
Sa {(2m + 1) • kgs) = iV • 4 [7r(2TO + 1)]"^. The structure 
factor of a uniform state takes a finite value at q — 
only: 5^(0) = N. This case can be regarded as the limit 
/i oo so that ko — and all the peaks of Sa{km) 
collapse into the peak at g = 0. 

B. Finite temperature 

Fig. [2] shows the two-point correlation functions ^ 
computed by MC simulations for a = 1.8 (domain ground 
state) and a = 3.0 (uniform ground state) at differ- 
ent temperatures (see Appendix [X] for details about the 
computational methods). In spite of the fact that the 
single-spin average ((cj)j)t is zero at any finite tem- 
perature, the correlation function reproduces the essen- 
tial aspects of the ground state spin configurations. For 
a = 1.8 (Fig. [2^), Ca{r) displays an oscillatory decay as 
a function of r, indicating that the loss of on-site mag- 
netization proceeds in such a way that the ground-state 
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FIG. 3: (Color online) Structure factor obtained by MC sim- 
ulations at T/g — 0.02 for a = 1.8 with Lorentzian fitting 
for the peaks located ai q — qm = (2m -|- 1) ■ qa,max (and 
J — 2.5). Inset: The HWHM Xa,m obtained by Lorentzian 
fitting of the MC data (full squares) and given by formula Q 
(line) is plotted versus q-m in a log-log scale, which clearly 
reveals the power-law behavior. 



segment-order is maintained. In the regime in which the 
ground state is uniform [a — 3), instead, the correla- 
tion function decays smoothly and, in general, monoton- 
ically (Fig. [2]d). A closer look at the highest reported 
temperatures {T/g = 0.8 , 1.4) reveals a small interval at 
short distances in which Ca {r) becomes negative (inset in 
Fig.Uja). This might be taken as an indication that, even 
starting from a uniform ground state, when the temper- 
ature is increased the system can spontaneously produce 
a phase with reduced symmetry in which the short-range 
order occurs with a well-defined modulation. We will 
come back to this point at the end of Section HITl 

In Fig. [HI the structure factor corresponding to a = 1.8 
and T = 0.02 is plotted. The set of discrete peaks of 
the ground state have broadened to Lorentzians centered 
at gr = {2m + 1) • qa,max- Here, qa.max means the po- 
sition of the highest-peak of the simulated Sa{q) at fi- 
nite T and does not, in general, coincide with kgs (the 
temperature dependence of qa,max will be discussed be- 
low). The occurrence of multiple peaks in the finite- 
temperature structure factor not only indicates that the 
periodic structure of the ground state propagates at fi- 
nite temperatures but also shows that some memory of 
the detailed square-wave spin profile is retained. As T is 
increased, peaks with m > rapidly lose weight and, for 
T > 0.1, basically only one peak is detectable. This im- 
plies a change of the correlation profile from triangular- 
wave-like (all harmonics) at low temperatures to cosine- 
like (single harmonic) at higher temperatures: the same 
cross-over is predicted to occur for the equilibrium mean- 
field spin profile within a 2d stripe-domain pattern and 
observed experimentally in the striped phase of ultra-thin 
Fe films grown epitaxially on Cu(OOl)^*'. Note that the 
height of the peaks of Sa{q) in the ground state scales like 
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FIG. 4: (Color online) Evolution of Sa{q) as a function of 
temperature: The red dotted line represents the peak position 
qa,max at different temperatures (parameters: J/g = 2.5 and 
a = 2.4, giving hgs — 51). Inset: Schematic view of a square- 
wave spin profile. 



(2m + 1)"^, while the ratio between the peaks at m = 1 
and m = in Fig. [3] is about one order of magnitude 
smaller at finite temperatures. In the next Section, we 
will give a simple explanation for this observation. The 
Lorentzian shape of the peaks and the m-dependence of 
their width (inset) - which are related to the exponential 
spatial dumping of the correlation shown in Fig. [5^ ~ will 
also be discussed in Section Hill 

A typical temperature dependence of the structure fac- 
tor is shown in Fig. [51 using a linear scale where only the 
most prominent peak m = is evident. Two facts are 
visible: 1.) the location of the maximum qa.max varies 
with temperature and 2.) the peak broadens consider- 
ably when the temperature is increased. We will discuss 
these two features more thoroughly. 



1. Temperature and a -dependence of qa.max 

(ha) = T^/qa.max IS plotted as a function of the tem- 
perature in Fig. [Slfor the set of values a = 1.6, 1.8, 2, 2.2, 
all in the regime a < for the chosen J/g. The 
ground-state value hg^^ found by minimizing the total 
energy (|B3p with respect to h, is also indicated. When a 
is increased - approaching the transition line to the uni- 
form state - both ground-state and finite-temperature 
values also increase. A strongly decaying long-range in- 
teraction favors longer periods. To be more quantitative, 
two temperature regions have to be considered: 

• For T/J >0.3, the period of modulation decreases 
with temperature, in a similar way to what is found 
for the stripe width in the Mean-Field Approxima- 
tion (MFA) of a similar but 2d model and in line 
with experimental results2£. 

• The temperature range T/J ^0.3 is more difficult 
to explain. The modulation period saturates at 
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FIG. 5: (Color online) Plot of (Jia) versus T/g in the domain- 
ground-state region. The computation parameters are L — 
1000, J/g — 2.5, and a = 1.6 (circles), 1.8 (squares), 2 (dia- 
monds), 2.2 (triangles). The statistical errors are smaller than 
the data symbols. 



the ground-state value for a — 1.6, 1.8 and remains 
below the ground-state value for a = 2.0, 2.2. 
We interpret the convergence of (/iq,) — > hgs with 
T ^ as a positive indication that our MC calcu- 
lations capture the essential equilibrium properties 
of the model, although we note that for larger peri- 
ods, in this temperature range, the MC acceptance 
rate approaches zero ("blocked condition"). A 
further investigation should be required to decide 
whether this is due to a technical limitation or 
rather to the set-in of intrinsic slow dynamics, by 
analogy with similar systemsS^i^i2^i2Si^. 

In Appendixiniwe will introduce an energy functional for 
finite T which depends parametrically on the period of 
modulation 2/i. Within some approximations, there we 
show that the minimum of such a functional is found for 
smaller and smaller h as the temperature is increased, 
thus reproducing qualitatively the dependence of qa,max 
on T. 



2. Temperature and a-dependence of the correlation length 
^a{T) = ^a\nax! ^a,max being the Half Width at Half 
Maximum (HWHM) of the Lorentzian centered at qa,max 



In Fig.m is plotted versus T/g for a = 1.6, 1.8, 2, 2.2 
(all falling in the region a < uq for J/g — 2.5) in a log- 
log scale. Dots correspond to MC data while the solid 
lines represent fits with the function Aa/T-^", with fit- 
ting parameters Aa and Ba ■ The best fit yields the same 
exponent Ba = 1.10(5) for each a, while Aa has a more 
complicated dependence on a, see squares in the inset of 
Fig. [SI (the zig-zag line will be discussed in the next Sec- 
tion). We conclude that the dependence of the correla- 
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FIG. 6: (Color online) Plot of versus T/g in the 
domain-ground-state region. The computation parameters 
are L =1000, J/g — 2.5, and a = 1.6 (circles), 1.8 (squares), 2 
(diamonds), 2.2 (triangles). The continuous lines represent 
the fit function y = Aa/x^" (see text). The statistical er- 
rors are smaller than the data symbols. Inset: The values 
of Aa/g obtained by fitting the results of the MC simulations 
(main frame) are reported as a function of a (full squares con- 
nected by dotted line). The theoretical value of the prefactor 
of l/T in formula ([11}, (4/7r^)/i'* {d^Egs/dh^), is indicated by 
the solid "zig-zag" line (see text). 



tion length on T is better described by Aa/T^°' than the 
Ising exponential relation e^/"^, which holds for g = 0. A 
deeper understanding of this difference will be provided 
in Section Hm For a comparison with the uniform regime 
(i.e., a = 2.6,2.8,3,3.2), let us consider the correlation 
function for a — 3.0 displayed in Fig. [^b. Note that the 
plot is limited to low enough temperatures in order to 
avoid the anomalous range of spatial decay where the cor- 
relation function becomes negative (see inset of Fig. [21d). 
Besides, in the temperature range T — 0.3... 0.8, the 
5 = behavior is recovered. In fact, looking at the cor- 
relation length ^Q, for a = 2.6,2.8,3,3.2, we find that it 
is better fitted by ~ exp(AQ./T), see the log-linear 
plot of versus g/T in Fig. [T] Remarkably, the energy 
barrier does depend on a, see squares in the inset of 
Fig. [7] (the dashed line will be discussed at the end of the 
next Section). 

Regarding the possibility to observe long-range order 
at finite temperature, our MC results seem to exclude 
such a hypothesis. In fact, in a correct analysis of the 
structure factor, beyond the usual intensive term con- 
nected with the correlations (Lorenzian-like functions), 
one should also take into account an extensive factor as- 
sociated with the occurrence of long-range order—. This 
last component has always been considered in the fitting 
procedure but it has never given a significant contribu- 
tion to the simulated structure factors ([3]). The occur- 
rence of long-range order at T = only is also supported 
by the low-temperature divergence of for both a < ao 
(Fig. [S]) and a > ao (Fig. [7]). The whole scenario con- 



FIG. 7: (Color online) Log-linear plot of versus g/T in 
the uniform-grourid-state region. The computation param- 
eters are L — 1000, J/g = 2.5, and a — 2.6 (circles), 2.8 
(squares), 3.0 (diamonds), 3.2 (triangles). The continuous 
lines are best fit functions. The statistical errors are smaller 
than the data symbols. Inset: The dependence of the barrier 
Ato = 2 J — 2g({a — 1) derived in Sect. Illll on a (dashed line) 
is compared with that of Ao, obtained by fitting computed 
in MC simulations (full squares). 



firms some recent theoretical works. In fact, in Ref. [T^ 
the absence of long-range order at any temperature for 
g > 0, J = and a > 1 is rigorously proved. Even 
if in that specific case a short-range ferromagnetic term 
was not included, it seems reasonable to extend such a 
result to the case J > and conclude that long-range 
order should not occur with the model ([1]) for a > 1—. 
The occurrence of a phase transition has been suggested, 
instead, for < a < 1 so that it would be particularly 
interesting to investigate the finite-temperature proper- 
ties of the model ([1]) in this regime. However, several 
other issues are related to the divergence of the energy 
per spin for pair-spin interaction decaying like 1 /r" with 
a < d, d being the dimension of the space in which the 
spin system is embedded, such as energy non-additivity 
and ensemble inequivalence^^^^ . For this reason our 
canonical MC method would not necessarily provide the 
unique and correct results in this context, but further 
analyses involving the comparison of different computa- 
tional approaches would be, most probably, required. 



III. DISCUSSION 

In this Section, we provide an explanation for some of 
the MC results on the basis of a simple physical model for 
the excited states of the Hamiltonian ([T]). In particular, 
we will provide a physical picture for the temperature 
dependence of the correlation length in the two distinct 
regimes a < and a> a^. 
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A. Case a < ao 

We construct excited states of the Hamiltonian ([T]) by 
modifying the square-wave profile to 



Sq(fco(j + 



am sin {km{j + Uj)) (5) 



with Uj being a displacement field. This perturbation 
corresponds to displacing the position of the wall be- 
tween adjacent segments, which creates a generally non- 
periodic spin configuration. The quantity we need to 
compute is the increment of energy due to the displace- 
ment field: 



A£h ^ £h[u] - £h[u = 0] 



(6) 



with £h being defined as {'H)j. A£h is computed 
perturbatively, i.e., in the limit of small Uq {uj 



setting fco 



q "9^ 



see Appendix [Cl For q « ko and 



kgs, with kgs being ir/hgs, one has 



A£a 



-E 



1 8'^£ 

L 2 '±±g± 2\~ |2 

o^gs 5, 2 9 1% 







(7) 



Eq. ([7]) describes the spectrum of the excited states 
(see also Refs. [ssUsgI for a model in 2d) and the coef- 
ficient of is a stiffness k'^^ {d^£gs/ dk^ against fluc- 
tuations from the ground-state spin configuration (see 
Appendix [C] for d^£gs/dk^ definition). Note the gapless, 
quasi-continuum nature of the spectrum of fluctuations, 
in clear contrast to the gapped spectrum of fluctuations 
in a pure {g = 0) Ising model. 

Eq. ([7]) is the central result of this Section, as it allows 
computing the structure factor Sa (q) and the correlation 
length (T) . The resulting structure factor ^ consists 
of a series of Lorentzian peaks centered at q = ±km 



1 °° f 

ni=0 I 



{q-kmf + Xl {q + kmf + Xl 



with a HWHM given by 



(8) 



(9) 



The reader is referred to Appendix [Cl for the details. The 
same behavior is observed in the MC results plotted in 
Fig. [3l Our analysis flnds the origin of the multiple peaks 
of Sa (q) in the quasi-continuum spectrum of gapless ex- 
citations (see Eq. ([7])) appearing in the frustrated model 
for a < aQ. A remarkable feature is the non-trivial scal- 
ing of the maxima with the higher-harmonic index 2m+\ 



Sa{q = ±km) 



1 ai 



2Xr, 



16 d^£gs 1 



7r2 5fc2 T{2m+iy 



, (10) 



which accounts for the strong reduction detected for the 
ratio between the peak heights for m — 1 and m — 
in the MC results at finite temperatures. Note also the 
square-power dependence of the HWHM Xa.,n on 2m + 1 
in formula ([9]). This theoretical prediction (solid line in 
the inset of Fig. [3] with log- log scale) is in excellent agree- 
ment with the behavior of the Sa{q) simulated for a = 1.8 
and J/g = 2.5 at T/g = 0.02 (squares in the inset of 
Fig. [3]). From the assumption ko = kgs (Eq. ^) it fol- 
lows that, within our analytic model, the highest peak of 
Sa{q) is expected to occur at qa^max — kgs and the cor- 
relation length is defined as = A~ g consistently. Even 
if we already know that in MC simulations qa,max does 
not remain constant as T is varied (see Fig. [5|) , this as- 
sumption produces a 1 /T-dependence of the correlation 
length 



1 



Xa,0 



d'^£ 

_ 2 9s 



1 



dkl 



T 



(11) 



which is in good agreement {B^ = 1.10(5)) with the cor- 
responding quantity computed again with the MC tech- 
nique (see Fig. [6|). Finally, the analytic model predicts 
that Aa — 2 (d^£gs/dkQ^ . Computing this expression nu- 
merically produces the solid curve in the inset of Fig. [6l 
The step-like behavior of 2 (d'^£gs/dkQ) reflects the fact 
that both the optimal domain width hgs and the second 
derivative of the energy ~ computed in /i = hgs - are dis- 
continuous functions of a'^ in virtue of the discreteness 
of the lattice. Both the order of magnitude and the scal- 
ing with a agree with MC calculations (squares in the 
inset of Fig. [5]) : Aa decreases as a increases approaching 
the uniform-ground-state region. In summary, the agree- 
ment between numerical and analytical results indicates 
that the distortion of the ground-state spin profile due 
to the displacement of domain walls represents the main 
disordering mechanism when a < a^. 
To the aim of reproducing the temperature dependence 
of Qa.max: tiiG expansion for q <^ kg — performed in Ap- 
pendix [C] to get from (|C12[) to (|C13[) - is not expected 
to be accurate anymore. However, in Appendix |D] we 
show that letting h be an adjustable parameter at finite 
T with an appropriate (temperature-dependent) stiffness 
we are able to reproduce qualitatively the decrease of 
the modulation period with increasing temperature ob- 
served in MC simulations. This, indeed, happens be- 
cause in the correlation function (jC21[) higher harmonics 
are progressively more suppressed as the temperature in- 
creases. As a result, the competition between the ferro- 
magnetic exchange and the antiferromagnetic long-range 
interaction turns out to be biassed with respect to the 
zero-temperature case and the period of modulation de- 
creases subsequently. This close relationship between the 
suppression of higher-harmonic components and the de- 
crease of the characteristic period of modulation has been 
already highlighted experimentally and by mean-field cal- 
culations in an equivalent 2d system, suggesting that it 
might be a general property of such models. 

The pure Ising Hamiltonian is invariant with respect to 
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any operation that changes the variable Uj to — (Tj-: it has 
the discrete symmetry group Z 2 . In the next Subsection, 
we will discuss this case in connection with a > a^. The 
l/T-dependence of the correlation length, obtained by 
introducing a long-range interaction {g ^ 0), suggests 
that, in the regime of a < ao, the frustrated system 
crosses over to the completely different universality class 
of one-dimensional chains hosting a planar spin field with 
SO (2) continuous symmetrjs^i^^. 

B. Case a > qo 

In this regime, the ground state is uniform and the 
Ising universality class is restored at low enough tem- 
peratures, as shown by the correlation length diverging 
exponentially as e"^/-^, see Fig. [T] Specific to this case is 
that A = Aq, see Inset Fig.[71 We try to explain this re- 
sult by considering that, in the pure Ising model = 0), 
the barrier A equals the energy cost to reverse half of 
the spins starting from a uniform configuration. Were 
the general arguments which associate such an energy 
with the low-temperature expansion of §^ applicable in 
the presence of long-range interaction, the energy of a 
single wall would be expected to equal A^. When half 
of the spins in the chain are reversed, the exchange en- 
ergy increases by 2J. To compute the variation due to 
the long-range interaction, note that this interaction en- 
ergy is just given by twice the interaction energy between 
the two parts of the chains lying on opposite sides with 
respect to the domain wall (as the self-energy in each do- 
main remains the same before and after the flip of half 
of the spins) . This interaction energy is given by 

A. = E = -2gE - -25C(« - 1) , 

(12) 

where C(a^) is the Riemann zeta function, while i and j are 
the site indices of spins lying on opposite sides of the do- 
main wall. The energy to create a wall becomes explicitly 
dependent on a and amounts to A^, = 2 J — 2gC,{a — 1). 
In the inset of Fig. [71 one can appreciate how this es- 
timate actually reproduces both the order of magnitude 
and the dependence on a of the energy barrier of the ex- 
ponentially diverging ^q, obtained from MC simulations. 
To be rigorous, one should point out that this approach 
is not completely justified in this context since, when a 
long-range interaction is present, the creation of a new 
domain wall is not statistically independent of the num- 
ber and the location of the pre-existing domain walls in 
the chain; such a hypothesis is indeed a basic assump- 
tion to put the correlation length in relationship with 
the cost to create a single wall in the system'^^ . Letting a 
go to infinity effectively reduces the spin-spin interaction 
to nearest neighbors only so that our system becomes 
equivalent to the usual Ising model, provided that the 
exchange interaction is replaced by J — g. 

At T = 0, the condition A„ = defines ao. In fact, as 



soon as A^, < the uniform configuration has no more 
the lowest energy and the system prefers to split into 
domains. For a given ratio J/5, ao fulfills the condition 
C(ao ~ 1) = J/ 9- Using the integral definition of the 
Riemann zeta function, the previous condition can be 
rewritten as 

— — C(ao — 1) = — ; / dx 

9 ° ^ r(ao-l)io e--l 



ao - 1 r(ao - 1) Jo (e" - 1)^ 

n^)Jo (l-e-2-)2' 

(13) 

this implicit equation for ao turns out to be exact^ (the 
solution being the solid line in Fig. [1]) . 

For completeness, we recall that at relatively high tem- 
peratures a well-defined period of modulation seems to 
emerge in the correlation function also for a > ao (see 
inset of Fig. [2)d). A naive, but essentially correct, inter- 
pretation of the temperature dependence ofqa,max in the 
regime a < ao suggests that thermal fluctuations effec- 
tively reduce the ratio J/g (the antiferromagnetic long- 
range interaction is fovored in the competition with the 
ferromagnetic exchange interaction, which finally leads 
to decrease the modulation period with respect to the 
T — case). In this sense, one may think that even 
when the uniform pattern has the minimum energy at 
T ^ (e.g. for J/g = 2.5 and a = 3 as in Fig. Hb), 
thermal fluctuations induce an effective decrease of the 
ratio J/g = 2.5 so that a modulated phase eventually has 
lower free energy at high enough temperatures. However, 
this effect can only be evident if such a crossover oc- 
curs when there is still enough correlation between spins 
to develop - at least ~ half-period of modulation, i.e. 
roughly for > l/qa,max- In fact, if the period of the 
underlying modulated phase is much larger than the cor- 
relation length ^a, two-point correlations just display a 
monotonic decay as a function of the lattice separation. 
A detailed investigation of this phenomenon would be, 
indeed, intriguing but it is beyond the purpose of the 
present work. 

IV. CONCLUSIONS 

The Mean-Field Approximation reported e.g. in 
Ref. [13 provides some straightforward results concern- 
ing ferromagnetic Ising system frustrated by a long- 
range interaction. However, the MFA fails in one im- 
portant instance: it predicts that the modulated order 
in the ground state propagates at finite temperatures up 
to a second-order transition temperature Tc, while the 
Landau-Peierls instability forbids a finite on-site {(Jj)T 
at any finite temperature^^i^. On the other side, MC 
simulations are much more accurate than the MFA, but 
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very difficult to perform under experimentally realistic 
conditions. For instance, the large modulation lengths 
often observed in experiments are practically inaccessi- 
ble to MC simulations. We concentrated on a model 
- Eq. Il]) - that is highly simplified but captures some 
essential characteristics of some physically relevant two- 
dimensional frustrated systems. Within this model, we 
have been able to enlarge the modulation length with re- 
spect to full two-dimensional MC simulations^ii^- . With 
this model, we have obtained a set of results that might 
help to shed light onto some experimental outcomes. 
In particular: The modulation length appearing in the 
ground state is found to remain a characteristic length 
at finite temperatures, where it appears as the length 
modulating the oscillatory part of the correlation func- 
tion. Strikingly, it decays with temperature in a way that 
is similar to the temperature dependence of the stripe- 
domain width observed in MFA and experimentally on 
Fe/Cu(001) films^*^i2i. In addition, the spatial profile of 
the correlation function contains the same kind of higher 
harmonics appearing in the MFA spin profile, with only 
one fundamental harmonic remaining at sufficiently high 
temperatures, as specified within the MFA and found 
experimentally^". In contrast to the MFA, which pre- 
dicts a second-order phase transition also in Id, we do 
not find any trace of a phase transition - and this is 
a major deviation from full two-dimensional MC simu- 
lations^ii^ or experimental findings. When the spatial 
decay of the long-range interaction is too short-ranged, 
the ground state and the finitc-tcmperature state lose the 
modulated character and become uniform. Correspond- 
ingly, the system crosses over from the universality class 
proper of Id systems with continuous symmetry^^i^ to 
the standard Id Ising-like universality class^i^. 
For future work, a more accurate treatment of the dis- 
placement field Uj beyond the q ^ kg approximation (see 
Appendix [C|) is certainly to be considered. 
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taining L spins and treat the long-range effects by repli- 
cating many identical copies of the "simulation box"iI. 
More explicitly, the interaction between two spins sepa- 
rated by r lattice sites reads 



G„(r) 



1 



1 



-nL[ 



(Al) 



where the index n accounts for the number N/ L of repli- 
cated boxes. Since we have in mind the thermodynamic 
limit N —^ oo, for numerical evaluation of Ga{r) we let n 
go to ±cx) in order to account for the copies of the system 
lying on both the left- and right-hand sides of the sim- 
ulated segment, containing just L spins. The effective 
coupling (|A1[) can be rewritten, in a way that is more 
suitable for computational purposes: 



G„(r) 




(A2) 

in the third passage we have neglected r/L with respect 
to M; the error of the whole approximation can be esti- 
mated following Ref. 113. This approximation reduces the 
main computational task to evaluating the finite sum over 
n, which is - however - rapidly convergent. Finally, the 
working Hamiltonian, restricted to our simulation box, is 
given by 



APPENDIX A: MONTE CARLO METHOD 

In this Appendix, we discuss the technical details of 
the MC method we used to study the finite-temperature 
properties of the Hamiltonian ([T]). A first important is- 
sue for the system under investigation is the treatment 
of finite-size effects. In the presence of long-range in- 
teractions, they need to be handled with particular care 
both numerically and analytically^^ . Some techniques to 
tackle the problem numerically are given, for instance, in 
Ref. l46l We perform our simulations on a system con- 



H = - J ^ a^a^+l + I E E ^'0-jGa(i - j) , (A3) 



.1 j=l 



which descends directly from ([T]) with the replica assump- 
tion ai±nL = <Ji {n = ±1 • ■ • ± oo), periodic boundary 
conditions on the simulation box ul+i = ai and setting 
r — \i — j\. Note that the indices i and j now vary in 
the range [1, L] and are allowed to be equal, Gq(0) being 
representative of the interaction between different spins 
in the original Hamiltonian ([T]); in this particular case 
{r — 0), there is no interaction inside the simulation box 
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but the z-th spin still interacts with its own copies lying 
in the different replicas, ai±nLi so that 



G„(0) 



^ 1 _ 2C(a) 

^ \nL\°' L°' 



(A4) 



The MC simulations have been performed using the 
Simulated Annealing (SA)^ paradigm. The SA is ex- 
tensively applied in statistical physics with the intent 
to study systems where both the ground-state energy 
and the equilibrium at low temperatures are inaccessi- 
ble through the basic Metropolis criterion^. Certainly, 
spin glasses^^, frustrated magnetic spin structures^ and 
models with long-range interactions^^ are some typical 
examples of systems where the SA and related methods^ 
are largely exploited. 

We also have to remind that in literature some clus- 
ter methods were employed in order to reach a correct 
thermodynamic equilibrium for a simple model where 
ferromagnetic long-range interactions are only presenlj^. 
However, the strong frustration due to the competition 
between the antiferromagnetic long-range interactions 
and the nearest-neighbor ferromagnetic exchange inter- 
action renders the generalization of such cluster MC tech- 
nique to the present case non-trivial. For these reason, 
we have followed in this work the main idea of Kirk- 
patrick et ali^. A random initial configuration (which 
should be considered as a paramagnetic state) is picked 
up. Subsequently, the thermodynamic equilibrium at a 
high enough temperature Tq is established. We remem- 
ber that the MC steps per spin considered here only com- 
prise Metropolis moves at the analyzed temperature. Tq 
is usually chosen in order to have a high MC acceptance 
ratio per spin. Then the temperature is decreased grad- 
ually T -> T - AT (AT > 0), and a fixed number of 
MC steps per spin r is run, starting with the last con- 
figuration sampled at the previous higher temperature. 
So, the main assumption is to force a constant and suffi- 
ciently slow cooling rate, defined as r = AT/t. We have 
taken AT = 0.1 . . . 0.001 and r = 1 . . . 5 x 10^ depending 
on the studied value of a. The procedure is completed 
when the ground state is approached. 

We have considered simulation boxes of size T=100, 
200, 500, 1000, and 2000. After discarding the first 
1x10^ MC steps, we have collected between 5x10^ and 
1 X 10^ measurements of the thermodynamic observables, 
repeating the simulation for each temperature at least 
three times. The estimation of the statistical errors 
has been achieved by applying the usual blocking tech- 

46 

mqu&i^. 



domains: a < ao. The lowest-energy configurations, at 
T = 0, are known to be given- by square- wave spin pro- 
files 

oo 

aj = Sq{koj) = ^ a„ sin {kmj) (Bl) 

m—O 

with ko = Tr/h, km = 7r(2m + l)/h and a„i = 
4/ [tt (2to 4- 1)]. With the orthogonality relation 
J2jLi e~^^^~^ = NSk^k' , the two-point correlations av- 
eraged over the site variables j can be computed: 



(Jj+r(Jj)j = J^YI ^'ii'^oU + r))Sq(koj) 

-j^ AT oo 

— ^ ^ a^' am sin {km{j + r)) sin {km, j) 

j — 1 m,m'— 

1 °° 

(B2) 

The Fourier coefficients of the series obtained in the final 
passage of equation (|B2p happen to be the same as for 
the symmetric triangular wave of period 2h so that in the 
text we use the compact notation {aj+r<^j)j = Tr(fco?'). 
Eq. (jB2[) allows writing the energy per spin for a general 
square- wave profile 



cos (kmr) 



-JTr(fco) + f E 4>E 



7n—0 r>l 



m=0 



(B3) 



which depends parametrically on the half-period of mod- 
ulation h. The ground-state energy for a given ra- 
tio J/g and a can be obtained by minimizing equa- 
tion (jB3p with respect to h numerically, which conse- 



quently defines the equilibrium domain width hgs at 
T = 0. The exchange term in (|B3p straightforwardly 
gives — ,/Tr(A;o) = ~J (1 ~ 2//i), also deducible by count- 
ing the number of walls present in the domain con- 
figuration with modulation period 2h. The function 
fa(km) = - {J/2)cos{km)+g/2Y,r>i (cos (fc^r) /r") in- 
troduced above will be used to write the perturbed en- 
ergy in a more compact form. 



APPENDIX B: CORRELATIONS IN THE 
GROUND STATE 



APPENDIX C: PERTURBATIVE TREATMENT 
OF CORRELATIONS AT FINITE 
TEMPERATURES 



In this Appendix, we compute the two-point correla- 
tions for a generic square-wave spin profile, representa- 
tive of the regime in which the ground state consists of 



In this Appendix, we develop a perturbative elastic 
model which allows us to compute the two-point corre- 
lations in the regime a < ao at finite temperatures. Let 
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us consider a displacement field, Uj, of the whole square- 
wave profile (|Bip : 



which can be exploited to average (|C3p with respect to 



aj = Sq {ko{j + Uj)) = ^ a,„ sin (fc„(j + uj)) . (CI) 



m=0 



To compute how the energy (|B2|) is modified by the pres- 
ence of this elementary perturbation, we introduce the 
constants 



a = krnU + r) 

b = k,n'j 



7 — k„iUj+r ^^^-j 
/3 — kjn'Uj 



where the two Greek letters will henceforth be assumed 
infinitesimal. Eq. (|B2p then involves terms like 

sin(a + 7) sin(6 + /3) = sin a sin b cos 7 cos (3 
+ sin a cos b cos 7 sin f3 + cos a sin 6 sin 7 cos (3 
+ cos a cos 5 sin 7 sin . 

(C3) 

We will further assume that the average over the lattice 
indices j, {. • . )j, can be performed independently for the 
rigid pattern variables (Latin letters) and for the fluctu- 
ating displacement field uj^: 

(sin a sin b cos 7 cos (3)j — (sin a sin b)j (cos 7 cos P)j . 

(C4) 

The average (. . . )j for elementary trigonometric func- 
tions with arguments a and b gives: 

{ (sin a sin 6) j = (cos a cos 6)^ = \6m,m' cos (fc^r) 
(sinacos&)j = — (cos a sin 6)^ = \5m,.m' sa\{kmr) , 

(C5) 



(sin(a -I- 7) sin(6 -I- /3))j 

= -^^m,m' COS (kmr) (cos 7 COS (3 + sin 7 sin /3)j 
+ -^Sm,m' sin (fcmr) (cos 7 sin /3 - sin7COS/3)j 



= 2'^™'™' {cos(fc„ir) (7?.e 
-|- sin (kmr) {Xm 



(C6) 



then, recalling that (3 — ^ — km [uj — Uj+r), we get: 



{aj+r'yj)j ^\Y1 {°^™ (cos(A:,„r) {He 



^ikm {Uj — 



sin (kmr) {Tm 



^ikm {Uj —Ujj^r) 



(C7) 



The introduction of the displacement field brings an in- 
crement to the energy of a general square-wave pro- 
file dm equal to 



^ OO 

A£h = -J-^{al^[sm{km)sm[k,n{uj-Uj+i)] + cos{kjn) {cos[kjn{uj - uj+i)] - 



ni—0 



m— r>l 

OO 



sin(fc™r) . r, , cos (/c„r) , r, , m 1^ 

Sm [k,n{Uj - Uj+r)\ H (cos [k,n[Uj ~ Uj+r)\ - 1) 



1 1 2 

J- ^ (a^i [- sin (km) kra{Uj+i - Uj) --COS (km) kl-^ [Uj+l - Uj) 



m=0 



m— r>l 



sin {kmr 



1 COS (A:,„r) 2 



^ km{uj+r Uj) 2 ^mWj+r Uj) 



(C8) 



where we have expanded the energy for small displace- To proceed in our calculation it is convenient to express 

ment differences Wj+r ~ the displacement field in terms of its Fourier transform 
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the averaged square difference is 



^3 iV ^ 



with 



, (C9) 



N 



Y,\u,\^\ - cos{qr)] , (CIO) 



the sum is performed over the Fourier wave numbers 
= ± (27rm) /N with m £ [-iV/2, 7V/2], but we drop 
the index m for simphcity. From Eq. (jC9|) it foUows that 



while (uj^ 



^01 3 



0. The previous results and 



the elementary trigonometric relation cos(a;) cos(y) 
(1/2) [cos(a; — y) + cos{x + y)] allow writing (jCSp as 



J 



cos(fc„i) - i cos(fc™ - g) - ^ cos(fc„ + q) 



\Un\' 



q 771—0 

~ E f E 1 "m^m^ E [cosik^r) -i cos [(fc„ - g)r] - ^ cos [(fc„ + g)r] 



— 



'2 ^ 



r 



(Cll) 



Recalling the definition of fa{km) (|B3p one can rewrite 
the perturbed energy (|C11|) as 



iv E E 

g m— 

1 



2 iL2 



2 /a (^n 



9) 



+ 2'^"(^™ + ?) ~ fa{km) 



for a 2d system with an analogous Hamiltonian. Within 
the range of validity of Eq. (jClSp and with restriction 
to h ~ hgs, an analytical formula for the structure fac- 
tor ^ can be derived. The thermal averages {■ ■ ■)t of 
the displacement field Uj, which appears in the perturbed 
two-point correlations (jC7[) . have to be performed first. 
Those thermal averages can easily be evaluated since the 
Hamiltonian (jClSp is quadratic for small perturbations of 
(C12) the ground state (h = hgs). The well-known theorem for 



As far as the large-distance behavior is concerned - like 
for the computation of the correlation length ~ one can 
expand the energy (jC12p for g <C fco to get 



iv E E 

q m—0 



^2 1 9'^ fa ^2\~ |2 



-9 \Uq 



(C13) 



where the derivatives are formally defined by assuming 
km to be a continuum variable k and taking the limit 
d^''fa/dk^ = hm/j^fe^ (9"/„/9A:"), which is, of course, 
more justified the larger h is. The fact that dkm/dko — 
km/ko and d^km/d^ko = implies 



9^ 

dki 



m=0 



m=0 



d'^ko dkm V / d'^k. 



krn 

ka 



dkl 



so that the perturbed energy (|C13[) finally reads 



-E 



]„1 ^^£h 2\~ |2 

2 dkl " ^' 



(C14) 



(C15) 



Eq. (jC15[) . specialized to ft. = hgs for the ground-state 
energy, essentially matches the result obtained in Ref . [1^ 



Gaussian distributed physical quantities^ readily gives: 



{{aj+raj)j)T = - ^ {a^^ cos (/c„ 



m=0 

X exp 



\k^{{{Uj -Uj+rf)j)T 



(C16) 



On top of the site average (|C10[) one has to perform the 
thermal average 

{{{u.+r ~ u,f),)T = ^ Y.(^1)t [1 - cos(gr)] ; (C17) 



in particular (wg)T can be computed ap- 
plying the equipartition theorem to (|C15[) : 
{kl/2) {d^£gs/dkl) {uDtQ^ = H. r/2 so that 



E' 



,2 N 



T 2 ^-^ 1 — cos(gr) 

,(C18) 

By writing the wave numbers explicitly, the sum in the 
previous formula can be evaluated analytically in the 
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thermodynamic limit (A'^ — > oo): 



z 1 ^ cos(qr) 2 



7V/2 



1 - COS (2^m) 



1 /27rr 



27r2 



TV 



m=-Af/2 

"^^^ 1 COS r^m) 

-^--^ 

771—0 m— 



TV 1 /27rr 



27r2 



N 



OC ^ OO 



COS 



TV 1 /27rr 



27r2 



iV 



m=0 

6 



TT / 27rr 

2 v~ 



1 /27rry 
iV 27r2r 



2tt^ N 



The thermal average (jC17p is finally obtained 

T 



{{{Uj+r - Mi) )j)T 



(C19) 



(C20) 



Combining (jC20p and (jC16[) . the sought-for quantity 
reads 



^ oo 

{{aj+r<Jj)j)T ^ T^Yl [flm COS (fc^r) e"-^°-'"'^] 



with 



T 



2-m 2 9^-£<,^ 



{2m + lf 



T 



(C21) 



(C22) 



The structure factor ^ is thus expected to have a series 
of Lorentzian peaks at g = ±(2to + 1) ■ fcgs, Xa^m being 
the corresponding HWHM. More explicitly 



m=0 I 



.m {q ~^ km) + 



(C23) 



APPENDIX D: OPTIMAL PERIOD OF 
MODULATION AT FINITE TEMPERATURES 

In this Appendix we provide a qualitative explanation 
for the dependence of qa,max on the temperature. First, 
we show that the decrease of the modulation period with 
increasing temperature is not reproduced just letting fcg 
be an adjustable parameter at any temperature. In fact, 
formula (jC21[) can be used to compute the two-point cor- 
relations associated with any square-wave profile, pro- 
vided that the appropriate stiffness against deviations 




FIG. 8; (Color online) Plot of {ha) versus T/g in the 
domain-ground-state region with J/g = 2.5 and a — 2: 
MC simulations (diamonds), elastic model with constant 
(dashed line) and temperature dependent (solid line) stiff- 
ness, kg (^d^ Eh/ d kg) and kg {'Hh)T /dk^l , respectively 

L J cos 

(see text). 



from the given period h ^ hgg is accordingly employed: 
fco [d'^Sh/dk^) (see (|B3P for the definition of £h)- In this 
way, one can account for the effect of thermal fluctua- 
tions on a square- wave profile of an arbitrary half-period 
h and construct the functional 



(Dl) 



where 



- CXD 

{{a,+r<Jj)j)T COS {k^r) e-^-"'-] (D2) 



with km = (2m + 1)tt / h {h =/= h„s are here allowed) and 
X^ m = kl^T/ [2fc2 (d^S^/dk^)]. The functional (Hh) 
(Eq. (|Dip ) can then be minimized with respect to h to 
obtain an effective equilibrium period of modulation at 
finite temperatures. This procedure produces the dashed 
line in Fig. [51 For J/g = 2.5 and a = 2, the optimal 
half-period of modulation corresponds the ground-state 
value, hgs — 11, for T < 0.07, while the functional (|Dip 
has a minimum in ft, = 10 for higher temperatures. All 
this indicates that the constant decrease of the modula- 
tion period observed in the MC simulations is not repro- 
duced just by including thermal fluctuations through a 
displacement field into the different square- wave profiles 
and by further minimizing the functional (|Dip with re- 
spect to h. Such a failure might be due to the assumption 
that the stiffness kg [d^Sh/dk^^ remains the same at any 
temperature. In order to circumvent this limitation, we 
propose a heuristic extension of our elastic model. Let us 
first compute the two-point correlations at an infinitesi- 
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mal temperature ST: 

{{a,+r<rj),)5T = ^ E K COS ikn.r) e-'^"-^] (D3) 



771—0 



with km = (2m + l)7r//i and SXa.m - 
kl^5T/ [2kl{d^£hldkl)]. The correlations dm 
can be thought of as resulting from a rigid spin profile 



CTj ^ a™ sin (g™j) 



(D4) 



777 = 



in which the wave numbers qm — (2m + l)qo are statisti- 
cally distributed. In particular, if a Lorentzian distribu- 
tion 



^^(9777) 



5Xa 



1 



(D5) 



is assumed, the corresponding averages - performed af- 
ter the site average (. . .)j - mimic the effect of thermal 
fluctuations such that Eq. (jD3p can then be rewritten as 



^ 00 



777=0 



dqmP{qm) COS {q 

(D6) 

The corresponding energy functional reads 



{'Hh)sT — {'Hh)q,„ = 



+ 00 



dqmP{qm) 



9 {^j+r^j)j 



E 



+00 



dqmP{qm)fa{q m ) 



(D7) 



To the aim of computing the correlation function at an 
infinitesimally higher temperature, the spin profile (ID4p 
can be further perturbed with a displacement field, which 
brings an increment to the energy functional (jD7p equal 
to 

(AHh),^ = E E "m / dqmP{qm) 



q m— 

^/q('?77i -q) + ^faiqm + q) - .fa{qm) 



(D8) 



By analogy with what done in the previous Section, we 
perform an expansion for q <^ kg (since q^s follow a 
Lorentzian distribution with maximum in fcoi Q ^ ^o 



as well): 

{AHh),^ - T7 E E "m / dqmP{qm) 
q ,77=0 



X 9m7: 



2 dql 



(D9) 

The fact that qm = (2m-|- 1)170 implies dqm/dqo = 9777/90, 
d^qm/d'^qo = and consequently 



-E- 

9 



dql 



(DIO) 



In the present case, the effective stiffness 
(qo {p'^'fi-h/dq^)q^ has a more complicated depen- 
dence on go with respect to Eq. (jClSp . However, we 
can simplify its computation significantly with the 
approximation 



c'9o 



2 /9' 



(DU) 

— c|,'^ meaning that the derivative with respect to 

(jKiQ I COS 

fco involves only the fluctuating functions, cos (kmr). The 
correlation function at the new temperature (T — ST + 
ST) is given by 



00 

777=0 



-t-00 



dqmP{qm) 



X cos (g^n,'') exp 



ST 



dk'^ I cos 



00 



+ 00 



dqmP{q ) cos (qmr) 



ST 



X exp 



2 1,2 dHH^) ST 

'^0 ai^ 



(D12) 

where in the last passage we have substituted qm inside 
the exponential with its maximum k'^. Such an approx- 
imation allows writing the energy functional at the new 
temperature again in the form (|D7[) . provided that the 
HWHM of the Lorentzian distribution P{qm) is changed 
into 



SXr, 



2 



ST 



ST 



JJ2 d^HH)5T I 

dk'^ I cos 



(D13) 



The whole process can then be iterated to obtain corre- 
lations at any temperature 



({(Tj+r(Tj)j)T 



1 °° 



cos (kmr) e 



-A„,„(T)r 



777=0 



(D14) 
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and the corresponding energy functional 

oo r ^ + 00 



ni—0 



with HWHM of P{qm) (letting ST dT) equal to 



(D15) 



2 r 

Jo 



1 fcj 

2 fc; 



2 (-T 



^0 



cJfcn I COS 



(D16) 

(Remember that the derivative with respect to kg in- 
volves only the fluctuating functions, cos(fcmr), and not 



the dumping terms). By minimizing numerically the 
functional (|D15p with respect to h we obtain the step- 
like curve in Fig. [H In this case a decrease with increas- 
ing temperature is indeed observed throughout the in- 
vestigated range. Such a qualitative agreement with MC 
results suggests that the change in the modulation pe- 
riod and in the effective stiffness, fcg d'^{Tih)T/dkQ , 

L J cos 

should be closely related. It is worth remarking that a 
better agreement is, probably, not to be expected given 
the expansion for g <C fco that we performed to pass 
from (jC12p to (|C13p and the further approximations in 
Eq. (|DTT|) and Eq. (|DT2|1 . 
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